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Negative X-ray reverberation time delays from 
MCG-6-30-15 and Mrk 766* 
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ABSTRACT 

We present an X-ray time lag analysis, as a function of Fourier frequency, for MCG- 
6-30-15 and Mrk 766 using long term XMM-Newton light curves in the 0.5-1.5 keV 
and the 2-4 keV energy bands, together with some physical modelling of the corre- 
sponding time lag spectra. Both the time lag spectra of MCG-6-30-15 and Mrk 766 
show negative values (i.e. soft band variations lag behind the corresponding hard band 
variations) at high frequencies, around 10 ~ 3 Hz, similar to those previously observed 
from 1H 0707-495. The remarkable morphological resemblance between the time lag 
spectra of MCG-6-30-15 and Mrk 766 indicate that the physical processes responsi- 
ble for the observed soft time delays is very similar in the two sources, favouring a 
reflection scenario from material situated very nearby to the central black hole. 
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1 INTRODUCTION 

The measurement of time lags as a function of Fourier 
frequency between soft and hard X-ray energy bands 
is widely used as a diagnostic of the emission mech- 
anism and source geometry in active galactic nuclei 
(AGN) (e.g. iPapadakis et alj l200ll; iMcHardy et"afl 120041 : 



.— ay • 

lArevalo et all l200fj . |2008l : ISriram et all 120091)" as well as 
X-ray binaries (XRBs) ( e.g. Miyamoto fc Kitamotd 1 19891 : 
iNowak fc Vaughan|[T99a : iNowak et all 1 19991 ). Usually they 
have positive values i.e. hard band variations lag soft band 
variations, and their origin is still ambiguous. Although posi- 
tive time lags are expect ed in standard comptonisation mod- 
els IjNowak et al.l 1 19991 ). they can also be caused by dif- 
fusive propagation of perturbations in the accretion flow 
(|Kotov et al.l(200ll ). 

Tentative detections of negative X-ray time lags 
(NXTL) (i.e. soft band variations lag hard band variations) 
were firstly reported fo r the type I S eyfert galaxies Mrk 76 6 
|Markowitz et al.ll2007l ) and Ark 564 (McHardy et al.ll2007l ). 
the latter commented that they could be due to reflection by 
material very close to the black hole ( BH). The only robust 
detection of NXTL, so far, came from iFabian et all (|2009l ) 
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for the narrow-line type I Seyfert galaxy 1H 0707-495. They 
detected NXTL of 30 s for frequencies greater than 6 x 10~ 4 
Hz. At lower frequencies, positive time lags, as commonly 
seen in AGN, were also detected. 

The orig i n of NXT L is, if anyth i ng, even les s clear. 
IFabian et alT (|2009h and IZoghbi et all (|2010l . l201ll ) favour 
the interpretation that they arise from reflection by matter 
lying very close to the BH. However, alternative explana- 
tions involving large scale distant reflectors with anisotropic 
geometry (e.g. accr etion disc wind) can also reproduce them 



(Miller et alThoiOl ). Thus, ideal NXTL candidates are the 
X-ray bright AGN which have long and continuous soft and 
hard X-ray light curves and which show evidence of repro- 
cessing e.g. possessing broad X-ray iron emission lines or 
having strong and variable soft excesses in their X-ray spec- 
tra. 

The type I Seyfert galaxy MCG-6-30-15 was the 
first AGN reported to possess a broad X-ray iron line 
ijTanaka et all fl995f) . Its reflection spectrum is considered 
archetypical for other AGN exhibiti ng similar broad iron 
line features (e.g. iNandra et afll2007l ). Based on the width 
of this line, using XMM-Newton observations, IWilms et aJj 
(2001) and IFabian et al.l (|2002l ) have put constraints for 
MCG-6-30-15 on the inner radius of the emission region 
of the order of 2 gravitational radii (r g ). Another X-ray 
bright type I Seyfert galaxy, observed repeatedly by XMM- 
Newton is Mrk 766. Although not one of the most convincing 
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cases, it shows weak evidence o f broad relativistic iron line 
(|Miller et al.ll2006l ; lMillerll2007f l. 

In this letter we re-analyse the archival XMM-Newton 
data of the above two sources and we report the existence 
of NXTL for both of them. Time lag analy ses of the same 
data sets have been already performed by IVaughan et al.l 
|2003l ) for MCG-6-30-15 and iMarkowitz et al.l (|2007h for 
Mrk766. These analyses used large time bins (100 s and 
60 s respectively) and were aimed to probe more the low 
frequencies part of the time lag spectrum (below 10~ 3 Hz). 
They measured positive time delays dropping reciprocally 
with the frequency. In addition both analyses used energy 
bands (0.2-0.7 vs. 0.7-2 and 0.2-0.5 vs. 0.5-1.1 respectively) 
not appropriate to enhance the existence of a possible reflec- 
tion signature. In Section [5] we describe the data sets and 
the data reduction procedures, in Section [3] we estimate and 
model the time lag spectra, and finally in Section [4] we dis- 
cuss our results. 



2 DATA SETS AND DATA REDUCTION 

MCG-6-30-15 has been observed six times and Mrk 766 eight 
times by XMM-Newton (all these data are publicly avail- 
able). To achieve the most accurate determination of the 
time lag spectrum in both sources, we selected the obser- 
vations which have an on-source exposure time larger than 
one full-day (i.e. more than 86.4 ks). The log for these ob- 
servations is listed in Table [T] 

The EPIC raw-data for all the observations are re- 
duced with t he XMM-Newton S cientific Analysis Sys- 
tem (SAS) (|Gabriel et all 120041 ) version 10.0.2. After re- 
processing the pn and the two MOS data sets with the 
epchain and the emchain SAS-tools respectively, we per- 
form a thorough check for pile-up using the task epatplot. 
For the production of the light curves we select events that 
are detected up to quadruple pixel-pattern on the CCDs i.e. 
PATTERN<12. For each observation and for each instru- 
ment separately, the corrected background-subtracted light 
curves in the 0.5-1.5 keV and 2-4 keV energy ranges are 
produced using the SAS-tool epiclccorr in time-bins of 20 
s. Then, for each observation and for each energy range we 
combine the background-subtracted light curves from pn, 
MOS 1 and MOS 2. Additionally, for each observation we 
also create the pn background light curve, in the 0.5-10 
keV energy range in order to trace time-periods with in- 
creased background activity that are going to be disregarded 
from our analysis. Typically, after the abovementioned data 
reduction procedures the length of the final effective light 
curves are of the order of 0.5-20 per cent less with respect 
to the initial pointing on-source time (see Table [T] columns 
'On time' and 'Eff.time'). 

Finally, for comparison, we also consider the 'at-least 
one-day long' publicly available archival XMM -Newton data 
for 1H 0707 -495. These are the ones used by iFabian et al.l 
<|2009l ) and IZoghbi et all (|201Cl . l201ll ) (see Table [1] for de- 
tails). 



3 TIME LAG ESTIMATION 

For each source and for each observati on we estimate, using 
the phases of the cross-spectrum fe.g |Prie"s tlcv 198lj), the 
time lag spectrum, r(/), at a given Fourier frequency, /, 
between the light curves in the soft, 0.5-1.5 keV, and hard, 
2-4 keV, energy bands. In order to optimise the contribution 
from the reflection component or the scattering/absorbing 
material, we choose the 0.5-1.5 keV band where the 'soft ex- 
cess', above an extrapolation from the higher energy power 
law continuum, is greatest. Below 0.5 keV the response of the 
PN instruments is not well calibrated. For the hard band, 
where we aim to isolate the continuum emission, we choose 
2-4 keV. This band is relatively free of contribution from the 
soft excess and also from the broad Fe Ka line. Moreover, 
above 4 keV Poisson noise becomes significant contaminat- 
ing further our data. During the estimation we average over 
a range of at least ten consecutive frequ ency bins and the 
corres pond ing error is compute d following iBendat fe Piersol 
l|l986h and lNowak et all l| 19991 ). The same procedure is re- 
peated for all the observations of the source, and we com- 
pute the average time lag value together with its error at 
the mean of the grouped frequency bins. 

The time lag spectrum of 1H 0707-495 (Fig. [[] right- 
ha nd panel red open-di a monds) has be e n alre a dy estimated 
by IFabian et all ll2009D; I Zoghbi et all (|2010l ); iMiller et al.l 
(|2010I ) and lZoghbi (■I a 1.1 [iOl lli and our result is in absolute 
accordance with them. Note that we use marginally different 
energy bands, with respect to the previous works cf. 0.3-1.0 
keV vs. 1.0-4.0 keV. 

In the left-hand panel of Fig. [1] we show the estimated 
time lag spectra for MCG-6-30-15 (filled black-circles) and 
Mrk 766 (open grey-diamonds) calculated using all of the 
data sets listed in Table [T] However, revolution 0999 of 
Mrk 766 has a much lower average count rate (i.e. much 
higher Poisson noise) than the other five observations, in- 
troducing considerable scatter to the time lag spectral esti- 
mates. Therefore, we recalculate the time lag spectrum for 
Mrk 766 without this observation and the results are shown 
in the same plot as filled red-diamonds. The time lag spec- 
tra in both cases are consistent within the errors, but after 
the exclusion of revolution 0999 the uncertainties are a good 
deal smaller. What is very interesting, is that the resultant 
time lag spectrum of Mrk 766 looks very similar to that of 
MCG-6-30-15. In both sources, the time lag spectra are pos- 
itive at frequencies below 6 x 10 -4 Hz, increasing towards 
lower frequencies with a power law like shape of slope around 
unity. At 6 x 10~ 4 Hz, the time lags in both sources become 
negative. They exhibit the same negative minima at around 
(7xl0 4 Hz, —18.9 ± 4.7 s) being inconsistent with zero at a 
99.99 per cent confidence level. They rise fast again to the 
positive time lag regime peaking at around (3 x 10 -3 Hz, 10 
s) and at even higher frequencies, they settle to a zero level, 
although a possible second excursion to the negative time 
lag regime (at around 10~~ 2 Hz) is possible. 

Given these similarities between MCG-6-30-15 and 
Mrk 766, after excluding revolution 0999, we combine the 
time lag spectra of MCG-6-30-15 and Mrk 766 to an aver- 
age time lag spectrum (ATLS hearafter). In the right-hand 
panel of Fig. [T] we show the ATLS together with that of 
1H 0707-495. In the same figure the two solid lines (black 
and red) correspond to the result of a least square fitting of a 
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Table 1. Summary of the XMM-Newton observations for all the sources. 



iJU HI t^L 


Obs. Id 


\-< 1 1 1 \~ i /~\n 

IXC V Ul Lit 


Olal I Llilll^ 

(UTC) 


I In i~ i yvi 1 1 

Vyil L1111L 

(ks) 


Rff time 
(k S £ 


VJUBCI Vlllfci 

mode^ 


VU LlllL laliC 




0029740101 


0301 


2001-08-01 16:07:44 


89.4 


81.3 


sw,sw,sw 


32.3±10.4 


MCG-6-30-15 


0029740701 


0302 


2001-08-02 04:25:17 


129.4 


123.3 


sw,sw,sw 


35.9±7.9 




0029740801 


0303 


2001-08-04 04:19:10 


130.5 


124.3 


sw,sw,sw 


32.9±13.5 




0109141301 


0265 


2001-05-20 08:33:34 


129.9 


113.6 


sw,fu,sw 


28.2±7.3 




0304030101 


0999 


2005-05-23 19:15:52 


95.5 


77.8 


sw,lw,lw 


5.8± 2.3 


Mrk 766 


0304030301 


1000 


2005-05-25 18:16:53 


98.9 


98.3 


sw.lw,lw 


11.8±4.6 


0304030401 


1001 


2005-05-27 18:19:03 


98.9 


93.1 


sw.lw,lw 


15.2±4.2 




0304030501 


1002 


2005-05-29 19:11:44 


95.5 


95.0 


sw.lw,lw 


18.6±4.5 




0304030601 


1003 


2005-05-31 18:12:53 


98.9 


98.4 


sw.lw,lw 


15.1±5.4 




0511580101 


1491 


2008-01-29 18:28:24 


123.8 


121.6 


lw,sw,sw 


3.0±1.3 


1H 0707-495 


0511580201 
0511580301 


1492 
1493 


2008-01-31 18:20:46 
2008-02-02 18:22:11 


123.7 
122.5 


102.1 
104.1 


lw,sw,sw 
lw,sw,sw 


4.6±1.6 
4.1±1.7 




0511580401 


1494 


2008-02-04 18:24:21 


121.9 


101.7 


lw,sw,sw 


3.2±1.4 



a The effective observing time, i.e. light curve duration, after the reduction of the raw data. 

^ For the pn, MOS 1, and MOS 2 respectively. sw:small-window, lwdarge window, and fw:fast-uncompressed. 

c The mean count rate and the standard deviation of the light curve in the 0.5-4 keV energy band. 




Frequency (Hz) Frequency (Hz) 

Figure 1. Time lag spectra for MCG-6-30-15, Mrk 766, and 1H 0707-495. Positive time lags indicate that variations in the 0.5—1.5 keV 
energy band lead the corresponding variations in the 2-4 keV energy band. (Left-hand panel). Time lag spectra for MCG-6-30-15 (black 
filled circles). Mrk 766 [excluding revolution 999] (red filled diamonds), [all the observations] (grey open diamonds). The dotted lines 
among the points of the time lag spectra are linear interpolations intended to guide only the eye. (Right-hand panel) Average time lag 
spectrum (ATLS) of Mrk 766 and MCG-6-30-15 (black filled circles) together with that of 1H 0707-495 (red open diamonds). The solid 
lines correspond to the best-fitting quadratic degree B-splines each one having ten control points. 



quadratic degree B-spline with ten control points to the two 
datasets, aiming only to outline the basic morphological dif- 
ferences/similarities between the two spectra. In comparison 
to the ATLS, the time lag spectrum of 1H 0707-495 seems 
to be shifted slightly to higher frequencies, remaining neg- 
ative at high frequencies and becoming broader around the 
deepest negative dip. The ATLS is negative over a smaller 
frequency range, and at higher frequencies exhibits a more 
pronounced 'oscillatory' behaviour. 



3.1 Reverberation and transfer functions 

In this section we attempt to fit the reverberation 
models that have been previously used for 1H 0707- 495 
ijFabian et al.ll2009l : lMnler et aLllioialZoghbi et alJl201lh to 
the ATLS. Since the uncertainties of the ATLS, at high fre- 
quencies (greater than 5 x 10 -3 Hz), are significantly smaller 
than those at lower frequencies, any best-fitting model will 
be driven by these high frequency data. Also above these 
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frequencies the Poisson noise becomes comparable to the in- 
trinsic source variability and the coherence goes significantly 
below unity (around 0.25) in which case the time lag mea- 
surements are no longer meaningful. Therefore, in order to 
quantify in a statistically robust way the physically inter- 
esting part of the observations i.e. the part of the time lag 
spectrum which shows the negative delays and the positive 
peak at around 3 x 10 -3 Hz, we fit the various models only 
to the ten lowest frequency points of our data set. Since the 
initial model-parameter values play a major role for this sort 
of nonlinear fitting their values are chosen in order to yield 
an initial model that is phenomenologically in accordance 
with the data. Note that the following best-fitting model re- 
sults are only indicative since none of the models possesses 
a single global minimum in the \ 2 hyper-surface. Also, the 
quoted errors are simply the diagonal elem ents of the covari- 
ance matrix (|Bevington fc Robinsonl ll992). 

In general, the time lag dependence on frequency can 
be used to constrain reflection models, and in particular 
the size of the refl ecto r. Following the pr evious studies of 
iMiller et all l|2010fl and lZoghbi et al.l (|201lh . we assume that 
a percentage / of the soft band flux is from a medium 
which reflects the continuum emission i.e. hard band flux, 
and that the reflector's response to the continuum is uni- 
form over a distribution of time delays, with a width of At, 
centered at to i.e. a top hat transfer function (Model 1). 
Such a model can approximate reflection from a spherical 
shell, from clumpy material not isotropically surrounding 
the source, and, to a first approximation, reflection from the 
disc around the central source. The maximum time delay in 
this case, t ma x = to + At/2, corresponds to the size of the 
reflector having a radius ci max /2. The \ 2 of the best-fitting 
model is 21.35 for 7 degrees of freedom (d.o.f.) and the best- 
fitting model parameters are At = 821 ±236 s, to = 887 ±75 
s and / = 0.13 ± 0.02. This is neither a statistically nor a 
morphologically acceptable fit, having a hypothesis proba- 
bility (n.h.p.) of 0.0033 and thus we do no show it. 

Then, following IMiller et all (|2010h . we considered the 
possibility that reverberation is likely to be present at dif- 
fering levels in both the soft and hard bands (Model 2) . We 
therefore considered two reflectors (i.e. two top hat trans- 
fer functions), of different size (i.e. different At), which 
physically corresponds to an extended reflecting medium, 
where the hard band photons travel further through the re- 
flecting medium, hence showing longer time delays. In this 
model the soft band transfer function is assumed to ex- 
tend from zero to t so ft (i.e. centered at t so ft /2) and the hard 
band transfer function from t so ft to t SOIt + Athard (i.e. cen- 
tered at t ao f t ± Athard /2). The \ 2 °f the best-fitting model 
is 13.18 for 6 d.o.f. (n.h.p=0.040) yielding t soft = 357 ± 113 
s, Athard = 1261 ± 105 s and two relatively small reflection 
fractions of / soft = 0.10 ± 0.04 and / har d = 0.15 ± 0.02. As 
we can see from the left-hand panel of Fig. O Model 2 seems 
to reproduce adequately the steep transition from the posi- 
tive to the negative time lag regimes as well as the negative 
minimum (points 4-7) but it does not reproduce well the 
subsequent positive peak (points 8-10). 

Finally, following IZoghbi et al.l (|201ll ). we consider a 
model which the hard band variations are caused by prop- 
agation of perturbations in the accretion flow and the soft 
lag variations are due to a refection component (Model 3). 
This model consists of a power law model for the positive 



hard lags and a top hat transfer-function for negative the 
soft lags. The \ 2 °f the best-fitting model is better, 9.63 
for 5 d.o.f. (n.h.p=0.0865) yielding for the soft component 
to = 142 ±45 s, At = 293 ±47 s and / = 0.23 ±0.10 and for 
the hard component a normalization of c = (4.4±1.8) x 10~ 3 
s and an index of q = —1.23 ± 0.18. As we can see from the 
right-hand panel of Fig. [2] Model 3 can very well reproduce 
the overall shape of the ATLS. 



4 SUMMARY AND DISCUSSION 

Using the longest XMM-Newton light curves in the 0.5-1.5 
keV and 2-4 keV energy bands we found that MCG-6-30-15 
and Mrk766 have remarkably similar time lag spectra and 
above 10 -3 Hz exhibit NXTL similar to t hose previously 
observed in 1H 0707-495 l|Fabian et al.ll2009l ,. 

With respect to our ATLS, if the NXTL arise through 
reverberation from scattering of radiation that is passing 
through partially opaque material (Model 2), the maxi- 
mum time delay of the order of 1600 s placing a limit on 
the extent of the reverberating region of approximately 800 
light-seconds (assuming a spherical geometry for the scatter- 
ing medium; this is the most reasonable assumption, even 
if the covering factor of the material is less than unity). 
This corresponds for MCG-6-30-15 (BH mass 5 x 10 6 M ; 
iMcHardv et all 120051) to 33 r g , and for Mrk 766 (BH mass 
1.8 x 10 6 M ; lBentz et al.ll2009h to 90 r s . Moreover, based 
on the best-fitting Model 2, the maximum delay for the 
soft band transfer function is 350 s corresponding to X- 
rays which are scattered while travelling through a medium, 
whose opacity decreases with increasing energy. Since, the 
soft band X-ray photons cannot travel far through the 
medium before being absorbed any scattered X-rays must 
have relatively short time delays. This value places an up- 
per limit on the region where the soft X-rays are scattered 
of 175 s, corresponding to a region of less than 7 r g and 20 
r g for MCG-6-30-15 and Mrk 766 respectively. Thus, even 
if the NXTL are due to scattering of the X-ray continuum 
from circumnuclear material, this material must be located 
within a region very close to the last stable orbit, at least in 
the case of MCG-6-30-15. 

We note that, the possibility of the continuum being 
scattered from clouds which are simply located close to the 
line of sight, but do not cover the sources isotropically, is 
probably ruled out by the vary fact that we observe NXTL 
in both MCG-6-30-15 and Mrk 766. Although the presence 
of such a specific geometry in one source, i.e. 1H 0707-495, 
could be justified as a chance coincidence, its presence in two 
additional sources renders this possibility rather unlikely. 

On the other hand, the high frequency shape of the 
ATLS (negative and then positive) rules out the possibility 
that there exists just a single time delay between the soft and 
hard X-ray bands, representing the light travel time between 
the X-ray source and the reflector. In this case we would 
expect NXTL at all high frequencies to have the same am- 
plitude, which is not the case in the ATLS (see Fig. [T] right- 
hand panel). The best-fitting Model 3, indicates a source 
size of the reflector no more than 150 s, corresponding to 6 
r g and 17 r g for MCG-6-30-15 and Mrk 766 respectively. 

In a top hat reverberation scenario, the wider frequency 
extent of the NXTL in 1H 0707-495 than in the ATLS im- 
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Figure 2. Best-fitting models for the ATLS. (Left-hand panel) Model 2 [Double top hat] (black thick solid line) together with the 
corresponding model components, soft (dotted line) and hard (dashed line). (Right-hand panel) Model 3 [Top hat with power law] (black 
thick solid line) together with the corresponding model components, soft (dotted line) and hard (dashed line). 



plies a smaller At and hence a smaller size of the reflection 
region. The smaller size implies a smaller BH mass, which is 
consistent with the overall shift of the lag spectrum to higher 
frequencies in 1H 0707-495. Conversely, the individual time 
lag spectra of Mrk 766 and MCG-6-30-15 are very similar; 
although the admittedly less precisely measured BH mass of 
MCG-6-30-15 is approximately twice that of Mrk 766. 

Regarding possible scalings to XRB s, the zero-crossin i 
frequency in GX 339-4 is around 4 Hz (|Uttlev et al.ll201ll 
which is only 5 x 10 3 times higher than those measured here. 
The maximum negative lag ( — 5 x 10" 3 s for GX 339-4 vs. 
-20 s from the ATLS) has a similar scaling, which is again 
less than the mass scaling, although reverberation may be 
less important in a hard state XRB. 

With respect to the power spectral density bend s 



7.6xl0" 5 Hz (MCG-6 -30-15; iMcHardy et all [20051 ). 



ly et 

4.4 xl0~ 4 Hz (Mrk 766: iMarkowitz et al l l2007f ) and 



1 .6 x 10"" Hz (1H 0707-495; IZoghbi et ap2010l ) they all seem 
to occur at lower frequencies where the time lags are strongly 
positive and they do not seem to scale in any simple manner 
with any particular frequency in the time lag spectra. 

Our analysis shows that NXTL are probably a common 
feature of AGN. However the simple models discussed here 
provide only an approximate description of the reprocessing 
mechanism. Therefore it is too early to speculate on the 
presence of any scaling law in the time lag spectral features. 



ACKNOWLEDGMENTS 

DE and IMM acknowledge the Science and Technol- 
ogy Facilities Council (STFC) for support under grant 
ST/G003084/1. IP acknowledges support by the EU FP7- 
REGPOT 206469 grant. Also, we would like to thank Phil 
Uttley for useful discussions. This research has made use of 
NASA's Astrophysics Data System Bibliographic Services. 



REFERENCES 

Arevalo P., McHardy I. M., Summons D. P., 2008, MNRAS, 
388, 211 

Arevalo P., Papadakis I. E., Uttley P., McHardy I. M., 

Brinkmann W., 2006, MNRAS, 372, 401 
Bendat J. S., Piersol A. G., 1986, Random Data: Analysis 

and Measurement Procedures. New York: John Wiley & 

Sons, — cl986, 2nd ed. 
Bentz M. C, et al, 2009, ApJ, 705, 199 
Bevington P. R., Robinson D. K., 1992, Data reduction 

and error analysis for the physical sciences. New York: 

McGraw-Hill, — cl992, 2nd ed. 
Fabian A. C, et al., 2002, MNRAS, 335, LI 
— , 2009, Nature, 459, 540 

Gabriel C, et al., 2004, in Astronomical Society of the Pa- 
cific Conference Series, Vol. 314, Astronomical Data Anal- 
ysis Software and Systems (ADASS) XIII, F. Ochsenbein, 
M. G. Allen, & D. Egret, ed., pp. 759-+ 

Kotov O., Churazov E., Gilfanov M., 2001, MNRAS, 327, 
799 

Markowitz A., Papadakis I., Arevalo P., Turner T. J., Miller 
L., Reeves J. N., 2007, ApJ, 656, 116 

McHardy I. M., Arevalo P., Uttley P., Papadakis I. E., 
Summons D. P., Brinkmann W., Page M. J., 2007, MN- 
RAS, 382, 985 

McHardy I. M., Gunn K. F., Uttley P., Goad M. R., 2005, 

MNRAS, 359, 1469 
McHardy I. M., Papadakis I. E., Uttley P., Page M. J., 

Mason K. O., 2004, MNRAS, 348, 783 
Miller J. M., 2007, ARA&A, 45, 441 

Miller L., Turner T. J., Reeves J. N., Braito V., 2010, MN- 
RAS, 408, 1928 

Miller L., Turner T. J., Reeves J. N., George I. M., Porquet 
D., Nandra K., Dovciak M., 2006, A&A, 453, L13 

Miyamoto S., Kitamoto S., 1989, Nature, 342, 773 

Nandra K., O'Neill P. M., George I. M., Reeves J. N., 2007, 



6 D. Emmanoulopoulos et al. 



MNRAS, 382, 194 
Nowak M. A., Vaughan B. A., 1996, MNRAS, 280, 227 
Nowak M. A., Vaughan B. A., Wilms J., Dove J. B., Begel- 

man M. C, 1999, ApJ, 510, 874 
Papadakis I. E., Nandra K., Kazanas D., 2001, ApJ, 554, 

L133 

Priestley M. B., 1981, Spectral Analysis and Time Series: 
Probability and Mathematical Statistics, Vol. 1-2. Lon- 
don: Academic Press, — cl981 

Sriram K., Agrawal V. K., Rao A. R., 2009, ApJ, 700, 1042 

Tanaka Y., et al., 1995, Nature, 375, 659 

Uttley P., Wilkinson T., Cassatella P., Wilms J., 
Pottschmidt K., Hanke M., Bock M., 2011, MNRAS, 
L250+ 

Vaughan S., Fabian A. C, Nandra K., 2003, MNRAS, 339, 
1237 

Wilms J., Reynolds C. S., Begelman M. C, Reeves J., 
Molendi S., Staubert R., Kendziorra E., 2001, MNRAS, 
328, L27 

Zoghbi A., Fabian A. C, Uttley P., Miniutti G., Gallo L. C, 
Reynolds C. S., Miller J. M., Ponti G., 2010, MNRAS, 401, 
2419 

Zoghbi A., Uttley P., Fabian A. C, 2011, MNRAS, 412, 59 

This paper has been typeset from a TgX/ DTgX file prepared 
by the author. 



